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State-by-state cross sections for dissociative electron attachment and electron-impact dissociation 
for molecular oxygen are computed using ab initio resonance curves calculated with the R-matrix 
method. When O2 is in its vibrational ground state, the main contribution for both processes comes 
from the 2 U U resonance state of but with a significant contribution from the 4 E“ resonant state. 
Vibrational excitation leads to an increased contribution from the low-lying 2 n g resonance, greatly 
increased cross sections for both processes and the threshold moving to lower energies. These results 
provide important input for models of 02-containing plasmas in non-equilibrium conditions. 
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I. INTRODUCTION 

Molecular oxygen is a major component of the Earth’s 
atmosphere and is fundamental for life. It plays an im¬ 
portant role in many natural and technological processes. 
In this paper we focus on the ability of low-energy elec¬ 
trons to dissociate molecular oxygen, in its fundamental 
electronic state, by means of resonant scattering. In par¬ 
ticular we consider the state-by-state processes of disso¬ 
ciative electron attachment (DEA), 

e" + 0 2 (X 3 E g -; v) -+ 0 2 "* -+ (T( 2 P) + 0( 3 P), (1) 

and electron impact dissociation (EID), 

e + 0 2 (X 3 Sg ; v ) —> 0 2 —t e +20 ( 3 P), (2) 

for each oxygen vibrational level v. The threshold for 
the process m is 3.64 eV, corresponding to the asymp¬ 
totic energy of O and O - fragments from the level v = 0, 
whereas for process © corresponds to the oxygen disso¬ 
ciation of 5.11 eV. 

Such collisions are important to study phenomena oc¬ 
curring in the upper atmosphere, re-entry physics, elec¬ 
trical discharges and plasma chemistry. DEA is the prin¬ 
cipal mechanisms, in molecular plasmas, for forming neg¬ 
ative ions from neutral molecules; the inverse process 
represents associative detachment. Although significant 
theoretical and experimental effort has been invested in 
characterizing electron-02 cross sections pj, information 
is only available for the vibrational ground state of 0 2 
for which both processes have rather small cross sections. 
Here, we will show that DEA and EID processes become 
much more important for vibrationally excited oxygen, 
as the corresponding cross sections increases by orders of 
magnitude with the excitation of the molecule. These 
cross sections represent fundamental input quantities 


’Electronic address: vincenzo.laporta@imip.cnr.it 


in kinetics models of oxygen-containing non-equilibrium 
plasmas, where high-lying vibrational levels of molecules 
can be hugely populated. 

The paper is organized as follow: In the Section HT1 we 
will give a brief account and numerical details on the 
adopted theoretical models for the calculation of the po¬ 
tential energy curves and on the description of the nuclear 
dynamics. The results are illustrated in Section Hill while 
the concluding remarks are given in Section IIVI 


II. THEORETICAL MODEL 

Many low-energy electron-molecule scattering pro¬ 
cesses are dominated by resonances. For molecular 
oxygen these processes can be quantitatively under¬ 
stood in term of four low-lying O^" resonances of sym¬ 
metry 2 n g , 2 n„, 4 X“ and 2 X“. First calculations 
were made by Noble et al. Q who used the R-matrix 
method [3] to determine the positions and widths for 
these four resonances. The 0 2 target was represented 
using nine states, corresponding to the orbital configura¬ 
tions: [core]l7r 4 l7r 2 and [core]l7r 3 l7r 3 . The scattering T- 
matrix was calculated, in the fixed-nuclei approximation, 
for internuclear distance range [1.85, 3.5] ao, using the 
configurations: [core]l7r 4 l7r 3 ( 2 n g ), [core] l7T 3 l7r 4 ( 2 n u ), 
[core]l7r 4 l7r 2 3cr u ( 4 S“) and [core]l7r 4 l7r 2 3cr„ ( 2 S“), for 
energies up to 15 eV. 

In our previous paper on state-to-state cross sections 
for resonant vibrational excitation (RVE) for electron- 
oxygen scattering Q, we extended Noble et alls cal¬ 
culations toward larger internuclear distances using the 
quantum chemistry code MOLPRO @. We used a multi¬ 
reference configuration interaction (MRCI) model and an 
aug-cc-pVQZ basis set. For the neutral ground-state of 
O2 as well as for the real part of its anionic state O^", 
the active space included the first 2 core orbitals (ler g , 
2<r g ) and 8 valence orbitals (2 a g , 3a g , 2 a u , 3 a u , I7 t u , 
l7r fl ) plus 150 external orbitals for a total of 160 contrac¬ 
tions. The MRCI orbitals were taken from a ground-state 
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multi-configuration self-consistent field (MCSCF) calcu¬ 
lation in which 4 electrons were kept frozen in the core 
orbitals and excitations among all the valence orbitals 
were considered for the 12- and 13-electrons of the neu¬ 
tral and ionic state respectively. 

In the present work we update Noble et a/.’s R-matrix 
calculations in order to have a better representation of 
the states involved in the scattering processes. Calcula¬ 
tions were performed using the UKRMol codes @. Or¬ 
bitals for the O 2 were generated using configuration inter¬ 
action (Cl) calculations using a cc-pVTZ Gaussian-type 
orbital (GTO) basis set. The electrons in the lowest 3 
core orbitals, (1 cr g , 2a g , ler u ) 6 , were frozen and an active 
space was constructed from 9 valence orbitals (3 a g , 2cr„, 
3a u , I 7 t u , 2 tt u , l7T g ) 10 . The scattering calculations with 
17 electrons also used a GTO basis [7| to represent the 
continuum electron. The calculations were based on use 
of a complete active space Cl representation which in¬ 
volves placing the extra electron in both a continuum or¬ 
bital and target orbitals using the following perscription: 
(1 <Jg, 2a g , lcr„) 6 (3cr g , 2 a u , 3 a u , l7r„, 2n u , l7r g ) n and 
(lcr g , 2cr g , lcr u ) 6 (3<r g , 2cr„, 3cr„, l7r u , 2ir u , l7r g ) 10 (4cr g , 
5er g , 37r„, 4a u , 2n g , 3n g , 1 <5 g , 1 d^) 1 , where the second set 
of configurations involves placing the extra electron in an 
uncontracted Q target virtual orbital. Calculations used 
an R-matrix radius of 10 ao- Both the MOLPRO and R- 
matrix calculations were performed using D 2 h symmetry. 

The potential curves calculated with the R-matrix and 
MOLPRO codes use different basis-sets, different theo¬ 
retical models and computational methods. This means 
that the resulting curves did not join and a merging pro¬ 
cedure was required. First of all the O 2 neutral ground- 
state potential was calculated using MOLPRO and all the 
subsequent results were referred to this curve. The real 
part of the resonance potentials, computed as a bound 
scattering state in the R-matrix calculations, were fixed 
at large internuclear distances beyond the crossing point 
to the MOLPRO curve which precisely reproduces the 
experimental oxygen atom electron affinity of 1.46 eV. 
At shorter internuclear distances, the R-matrix resonant 
curves, are simply a smooth continuation of these large 
R curves and therefore their positions are fixed by the 
corresponding MOLPRO curve. Potential energy curves 
for the 02(X 3 S g ) ground state, which supports 42 vi¬ 
brational states, and the four resonance states, plus the 
corresponding widths, are reported in Fig. [L] We note 
that the new calculated potentials are slightly changed 
with respect to those of our previous paper 0, in partic¬ 
ular for the 2 n„ symmetry, and this could affect the RVE 
cross section values. To check this aspect, we performed 
new calculations for RVE processes, observing however 
no significant changes in the cross sections with respect 
to the corresponding results in 0. 

Previously 0 we showed that for energies up to about 
4 eV the RVE cross sections are characterized by narrow 
spikes dominated by the 2 II g resonance. For energies 
about 10 eV, the cross sections showed a broad maxi¬ 
mum due to enhancement by the 4 X“ resonance with 
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FIG. 1: (Color online) (a) Potential energy curves for ground 
state X 3 S~ of O 2 and for the four lowest resonant states O^"; 
(b) The corresponding resonance widths, r (R), as a function 
of the internuclear distance. 

a smaller contribution coming from the 2 E“ state. In 
the case of DEA and EID processes we show that they 
are dominated, at least for low vibrational levels v of 
O 2 , by the 2 H U resonance. Although electron attach¬ 
ment from oxygen molecules has been widely studied 0 
there are only a few, rather old, papers reporting cross 
section measurements of the DEA process as a function 
of the incident electron energy. We cite here Rapp and 
Briglia 0 , Schulz 0 and Christophorou et al. fl2l . 
DEA from the 1 A g state of O 2 has also been observed jl3[ 
and found to have a cross section of similar magnitude to 
process Q. Cross sections for EID only appear to have 
been measured for electronically excited states of O 2 by 
Cosby 0. 

In Born-Oppenheimer approximation, the dynamics 
of the oxygen nuclei in DEA is treated within the 
local-complex-potential model 0, already satisfactory 
adopted for resonant calculations in other diatomic 
molecules (see Refs. 0,0 and references therein). The 
corresponding cross section from the vibrational level v 
of oxygen and for electron energy e is given by: 

cr v (e) = 2n 2 ^ — lim |£(I?)| 2 , (3) 

K /1 R-¥ 00 

where K is the asymptotic momentum of the dissociating 
fragments O and O - with reduced mass /z; m e and k = 
y/2m e e are the incoming electron mass and momentum 
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respectively and £(i?) is the solution of the Schrodinger- 
like equation for the resonant state and total energy E = 
e v + e: 


where V~ + is the complex potential of the resonance 
reported in Fig. |U V 2 = T/(2nk) is the discrete-to- 
continuum potential coupling and Xv is the vibrational 
wave function of O 2 corresponding to the vibrational level 
v. R represents the internuclear distance. There is not in¬ 
terference among the resonant states with different sym¬ 
metries and for the two £“ states the interference can 
be considered negligible as they have different spin mul¬ 
tiplicity. 


III. RESULTS AND DISCUSSION 



Electron energy (eV) 


Figure^a) shows our results for the DEA cross section 
for O 2 (u = 0) compared with the experimental results of 
Rapp and Briglia [10], Schulz El and Christophorou et 
al. [121 ]. The agreement is excellent, within the experi¬ 
mental error of ±15%: the peak is positioned at 6.43 eV 
with an absolute value of 0.0154 A 2 and the FWHM 
of about 2 eV. Fig. HDb) reports cross sections coming 
from all four resonances and their sum. As found exper¬ 
imentally, the main contribution to DEA cross section 
comes from the 2 H U state of 0^~. We also note, at about 
8.5 eV, the presence of a significant contribution due to 
the 4 £“ symmetry which, added to the main 2 II„ contri¬ 
bution, reproduces the high-energy tail observed exper¬ 
imentally. The non-negligible contribution of 4 £“ res¬ 
onance to DEA cross section has already been deduced 
from the measured angular distribution of O” ions [l6j ] . 

Figure [3] shows our DEA cross sections calculated for 
higher vibrational levels of O 2 , v = 10,20 and 30, com¬ 
pared with the result v = 0. All cross sections are 
summed over all four resonance states. As expected 53, 
the threshold of the process shifts to lower energies as 
the vibrational level increases; the is due to the reduced 
threshold for dissociation limit. At same time the maxi¬ 
mum value of the cross section for v < 30 grows by of sev¬ 
eral orders of magnitude at low en erg ies. This behaviour 
is due to the survival factor e~ p [l8j ]. given approxima¬ 
tive^ by: 


e p ~ exp 


f Rr T(R)dR 
Re hv{R) 


(5) 


where integration is extended over the region between the 
classical turning point ( R e ) and the stabilization point 
(Rc) and v(R) is the classical velocity of dissociation. 
As the vibrational level increases, R e increases and the 
survival probability grows like the cross section. The os¬ 
cillations, visible at high vibrational levels, results from 


FIG. 2: (Color online) (a) Calculated dissociative electron 
attachment cross section for v = 0 compared with exper¬ 
imental results of Rapp and Briglia 0 , Schulz 0 and 
Christophorou et al. [12|. (b) Contributions from the four 
resonances to the total cross section. 


the interplay between the neutral vibrational wave func¬ 
tion and the resonant continuum wave function. Anal¬ 
ysis of the resonance contributions to the cross section 
at v = 20 shows a complicated picture with both II res¬ 
onances making significant contributions at low energy 
and important contributions from the £ resonances at 
higher energies. 

Figure [3] also shows the cross sections for v = 30. The 
increasing trend of the maxima is now inverted. The 
eigenvalue for this vibrational level Q, in fact, lies above 
the energy of the 0“( 2 P) ± 0( 3 P) asymptotic state (see 
Fig. so that the DEA process becomes exothermic 
with a threshold at zero-energy. For this case R c > R e 
so that the exponent p in Eq. © vanishes (T(R > R c ) = 
0) and the magnitude of the cross section is no longer 
governed by the survival factor. The decreasing trend of 
the cross section maxima is maintained for v > 30. 

The corresponding EID cross section from 
C> 2 (X 3 £“;?;) and electron energy e is given by fl9| : 

^ EID W = J de'j \( X AR)\V\m)\ 2 , (6) 

where (• • •) means integration over the internuclear dis¬ 
tance R , £(i?) is the resonant wave function solution of 
the Eq. (HD and xv is the continuum wave function of O 2 
with energy e! representing the 20 + e~ fragments. The 
continuum energy e' was integrated from O 2 dissociation 
threshold up to 10 eV. 

Figure 0Da) shows the calculated EID cross sections 
for v = 0. Contributions coming from the four O^" res¬ 
onances are shown. Like DEA, the main contribution 
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FIG. 3: (Color online) Calculated dissociative electron at¬ 
tachment cross sections for vibrationally excited O 2 , for lev¬ 
els v = 10, 20 and 30, compared with the ground level v = 0. 
The cross sections include all the four resonant contributions. 


comes from the 2 II U symmetry, which gives a maximum 
at 6.93 eV with an absolute value of 0.0072 A 2 , and a sig¬ 
nificant contribution from the 4 S“ at 9.47 eV with a cross 
section of 0.0026 A 2 . The contribution from the other 
symmetries is negligible. Figure IHb) shows the cross 
sections for excited vibrational levels of oxygen compared 
with those for the ground level. Cross sections include all 
four resonance contributions. As expected the EID 
cross section increases rapidly with vibrational excitation 
and the threshold moves to lower energy. The cross sec¬ 
tion for v = 30, for example, reaches a value comparable 
with that of the well-known Shuman-Runge dissociative 
transitions for the same vibrational level [20( , whose dis¬ 
sociative cross sections decrease with increasing the vi¬ 
brational excitation of the molecule. We note that EID is 
known to occur efficiently via electron impact excitation 
of o 2 [Ml. However the threshold for electronic exci¬ 
tation processes lie significantly higher than those con¬ 
sidered here, making them relatively unimportant in low 
temperature plasmas. 


IV. CONCLUSIONS 

In conclusion, we present new calculations for dissocia¬ 
tive electron attachment for oxygen using ab initio po¬ 
tential energy curves and the first calculations for reso¬ 
nant electron-impact dissociation of oxygen. We confirm 
the dominant contribution of 2 H U symmetry in both pro¬ 
cesses starting from v = 0. Both cross sections however 
increase rapidly with increase in the vibrational excita¬ 
tion of the molecule and for these vibrationally excited 
states it is necessary to consider the contributions from 
all four of the low-lying O^" resonance states. 

Finally, Fig. [5] compares cross sections for different 
resonant processes involving electron-0 2 (X 3 £“, u) scat¬ 
tering: ourprevious results for v —> v + 1 vibrational- 
excitation Q, with the present cross sections for DEA 
and EID for the cases when the 0 2 is initially in its vibra- 
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FIG. 4: (Color online) (a) Calculated cross sections for reso¬ 
nant electron impact dissociation of oxygen for v = 0. Con¬ 
tributions from the four O^" resonances are shown; the black 
curve is the sum. (b) Dissociating cross sections for vibra¬ 
tionally excited oxygen at v = 10, 20 and 30 compared with 
the results for v = 0. 


tional ground state and in its v = 20 state. For the v = 0 
level vibrational excitation is the dominant low-energy 
process. However for vibrationally excited molecules the 
DEA increases rapidly and becomes much the most im¬ 
portant process at low energies. We expect that the in¬ 
clusion of these results will have important consequences 
in models of plasmas containing molecular oxygen and, 
in particular, will lead to a significant increase in O - ion 
production in these models. 

Recently, we have calculated also the rate coefficients 
for the processes m and © as a function of the electronic 
temperature and for all the vibrational levels. Details of 
the calculations and results will be reported in [2lJ. DEA 
and EID cross sections as a function of electron energy 
and 0 2 vibrational state, along with the corresponding 
rate coefficients, have been made freely available in the 
Phys4Entry database [22j . 
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FIG. 5: (Color online) Summary of the cross sections for 
electron-02 collisions: Resonant vibrational excitation for 
Av = 1 transition [j (solid line), dissociative electron at¬ 
tachment (dashed line) and resonant dissociation by electron 
impact (dotted line) for (a) v = 0 and (b) v = 20. The cross 
sections are the sum over the four Of resonances. 
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